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. The compatibility of neutrino-nucleus deep inelastic scattering data within the universal, 

factorizable nuclear parton distribution functions has been studied independently by several 
^ ' groups in the last few years. The conclusions are contradictory, ranging from a violation 



of the universality up to a good agreement, most of the controversy originating from the 
use of the neutrino-nucleus data from the NuTeV collaboration. Here, we pay attention 
to non-negligible differences in the absolute normalization between different neutrino data 
sets. We find that such variations are large enough to prevent a tensionless fit to all data 
simultaneously and could therefore misleadingly point towards non-universal nuclear effects. 
We propose a concrete method to deal with the absolute normalization and show that an 
agreement between independent neutrino data sets is established. 

o, 

The global analyses of the nuclear parton distribution functions (nPDFs) study the applicability 

(N ■ 

| of the collinear factorization in hard processes involving bound nucleons. The most recent analyses 

^ . [l]-4] include data on charged lepton nuclear deep inelastic scattering (DIS); Drell-Yan dilepton 

production in proton-nucleus collisions; and, in some cases, production of high transverse momen- 



tum pions in deuteron-gold collisions and neutrino-nucleus DIS. The good overall agreement with 



the available high energy data supports the existence of universal, process independent nPDFs. 

The neutrino- nucleus DIS, frequently used also in constraining the PDFs of the un-bound pro- 
ton, has recently invoked a special attention as its compatibility within the framework of universal 
nPDFs was questioned 0] . It was even declared [f| that there is no way to satisfactorily repro- 
duce the neutrino- nucleus and the other nuclear data simultaneously with a single set of nPDFs. 
This could have far-reaching consequences as the inability to to find a set of nPDFs which at the 
same time describes all the considered data is the expected sign e.g. of a violation of the univer- 
sality, or a breakdown of the DGLAP evolution. However, the same signal can also occur if one or 
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more of the experimental data sets contain unrecognized systematic inaccuracies. 



Contradictory results were first presented in 



where up-to-date nPDFs were found to give 



an excellent overall agreement with neutrino data from CDflSW j|, CHORUS [l(J and NuTeV j3] 
collaborations, although issues with the normalization of the NuTeV data were identified possibly 
explaining the results of [a]. Similar conclusions were reached in J] where data from all these 
experiments were utilized in a global nPDF analysis without an apparent disagreement. However, 
the baseline PDFs utilized there [ill were already constrained by the NuTeV D1S data and their 
uncertainties were treated as additional, uncorrelated point-to-point errors. Furthermore, the anal- 
ysis did not use the absolute cross-sections, but the far more scarce structure function data. Given 
all this, the neutrino data did not carry as heavy importance as in {(J. For more comprehensive 
review of the present situation, see [l^ . 

In the present article, we will show that when accounting for the overall normalization of 
the experimental data in neutrino DIS, all three data sets do show a uniform pattern of nuclear 
modifications, well reproduced by the existing nPDFs. This reinforces the conclusions of [§], in 
a model-independent way, supporting the functionality of the factorization in neutrino DIS. We 
make the point even more concrete by employing a method based on the Hessian error analysis to 
verify the consistency of these data with CTEQ6.6 Q and EPS09 [1] global fits. 
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FIG. 1: The neutrino and antineutrino data presented as -^Average Q e ^> panels), and as -^Average (right 
panels). The CHORUS (blue circles) and CDHSW (green diamonds) data has been horizontally shifted 
from the NuTeV (black squares) data points for clarity. 



We utilize the neutrino-nucleus DIS data from the NuTeV [J, CHORUS Q and CDHSW 
Q] experiments. The difficulty in dealing with the neutrino data is that no reference data from 
hydrogen or deuterium target are available and we are forced to use the absolute experimental 
cross-sections cjg xp instead of cross-section ratios. However, in order to better see the nuclear 
effects we still prefer to present the data as ratios 

&*™(x, v, E) 

lT(x,y,E)= "* { g J (1) 

°CTEQ6.6^ X ' f: ^ > 

where the theoretical cross-sections Cq TE q 6 6 are calculated with the CTEQ6.6M central set. As in 
the theoretical calculations include corrections for the target mass and electroweak radiation, 
and are carried out in the SACOT-prescription [14] of the variable flavour number scheme. In 
order to avoid higher-twist effects we restrict the virtuality Q 2 and the final state invariant mass 
W by conditions Q 2 cnt > 4GeV 2 , and W 2 ut > 12.25 GeV 2 . This leaves us with 2136 NuTeV, 824 
CHORUS, and 937 CDHSW data points. For a concise presentation of this large amount of data, 
we form an average 

\«6hxcda; / \igfixedx / \j6fixcdx / 

where <5j is the experimental error (statistical and systematic added in quadrature), and the sum 
runs over all data points in the same x-bin. This procedure neatly summarizes the main features 
of the neutrino data as a function of x but we stress that it is used here only for plotting the data, 
the numerical results being computed using the absolute cross-sections. The ratios constructed 
this way are shown in the left panels of Figure [TJ Although the data from different experiments 
appear to be in rough mutual agreement, the scatter is still non-negligible. Especially the NuTeV 
neutrino data seems to lie systematically below the rest and as such is likely to trigger tension in 
a global fit — especially so if the NuTeV correlated systematic errors are taken seriously as in 



171 ] . However, as a function of x the shape of the data seems to follow the usual nuclear effect 



suggesting that the problem is rather in the absolute normalization, as already conjectured in 
For this reason, we define 



8]. 



/ cxp( £: ) = a expA X iVi E ) X B i( x ,V), ( 3 ) 
iG fixed E 

and similarly for the theoretical calculation. The factor Bi(x,y) represents the size of the ex- 
perimental (x,y)-bin making I£ (E) thereby an estimate for the integrated cross-section. Now, 
instead of Eq. ([T]) we consider the ratio of the normalized cross-sections 

R ( *' y ' E) = TP — (t 7, f\it» — m- (4) 

°CTEQ6.6^ X ' yi ^ll 1 CTECl6.6\ rj ) 
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The averaged neutrino and antineutrino data normalized in this way are plotted in the right-hand 
panels of Figure [H demonstrating how all the considered data seem to fall in agreement. Espe- 
cially, the NuTeV neutrino data has moved upwards while the CHORUS and CDHSW neutrino 
data have remained essentially unchanged. This observation suggests that the origin of the diffi- 
culties in accommodating the neutrino data in a global fit [(| is due to an unnoticed problem in 
the experimental normalization of the NuTeV data — that the uncertainties have probably been 
underestimated by the experiment. 
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FIG. 2: The experimental i?Avora g c compared to the predictions from CTEQ6.6 and EPS09. 

In order to see how the normalized data compares with the predictions from the present nuclear 
PDFs, we replace in Eq. (J3]) the experimental cross-sections by the theoretical ones computed with 
the bound proton PDFs ff~(x,Q 2 ) obtained standardly by 

/^(x,Q 2 ) = ^ EPS09 (x,g 2 )/f TE Q 6 - 6M (,,Q 2 ), (5) 

where the factor ^^' EPS09 represents the EPS09 [lj] nuclear modification in free proton PDF 
/ CTEQ6.6M (a . jQ2) . The results are shown in Figure [2j where the data points are the same as 
in the right-hand panels of Figure [H and the band indicates the theoretical calculations with all 
PDF uncertainties added in quadrature Q|. The good agreement indicates that it should be pos- 
sible to include these data in global fits without significant mutual disagreement or tension with 
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the other data sets. We note that in the normalization procedure described here, also part of the 
PDF uncertainties cancel thereby making the theoretical predictions more solid. 

We turn now to a more quantitative description of the data sets accounting for the normalization. 
The technique described here is based on the Hessian uncertainty analysis 



151 ] performed e.g. in 

the EPS09 and CTEQ6.6 global fits [la ]. The neighborhood of the minimum x 2 1S approximated 
by an expansion 

X 2 « xl + Sa i R ij Sa j = X 2 o + ^2z 2 , (6) 

ij i 

where 5dj is the deviation of the fit parameter a,j from its best-fit value. By diagonalizing the 
Hessian matrix Hij one finds the uncorrelated parameter directions Zj in terms of which the central 
set So, and the error sets are defined: 

z(S ) = (0,0,..,0) 



z(Sf) = ± V / Ax2(l,0,...,0) (7) 



z(St) = ±VA^(0,1,...,0) 



where A\ 2 is the maximum permitted deviation from the minimum x 2 ■ These sets enable the 
calculation of any PDF-dependent quantity X at the origin and at the corners of the z space, but 
in order to obtain an estimate in an arbitrary point z(S) = (zi, Z2, ■ ■ ■) close to the origin, we need 
to use a linear approximation 

x[s]^x[s ] + Y, dX[S] 

k 

where 



dz k 



s=s 



Zk 



X Q + D • w, 



(8) 



x [st] - x [s k 



Zk 



(9) 
(10) 



Let us now consider a larger data set {X data }. The agreement with the PDF set S can be quantified 
by formally adding its \ 2 contribution x 2 {^ data } to Eq. ([6]) 



x 2 = xl+ Yl 



X k [S]-X^' 2 



rdata 
°k 



(11) 
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NuTeV 
Normalization 
No normalization 

CHORUS 
Normalization 
No normalization 

CDHSW 
Normalization 
No normalization 



All CTEQ6.6 and EPS09 error sets 



xl=o/N Xw miJl /N EPS09-penalty CTEQ-penalty 



0.84 
1.04 



0.77 
0.90 



13.9 
40.3 



35.4 
42.5 



xi=o/ N xlnin/N EPS09-penalty CTEQ-penalty 
0.70 0.69 2.13 2.63 

0.86 0.81 3.35 14.4 

xl=o/N Xw miD /N EPS09-penalty CTEQ-penalty 
0.70 0.64 7.20 17.3 

0.81 0.74 10.4 17.8 



Only EPS09 error sets 



airi /N EPS09-penalty 



0.81 
0.94 



33.8 
77.4 



Xw ■ JN EPS09-penalty 



0.70 
0.84 



2.48 
5.13 



Xi min /N EPS09-penalty 



0.68 
0.78 



9.26 
14.1 



TABLE I: The \ 2 /N values for the neutrino data and the penalties induced in EPS09 and CTEQ6.6. The 
left block of the table corresponds to the analysis including both EPS09 and CTEQ6.6 error sets, while the 
right block corresponds to keeping the free proton PDFs fixed at their central value. 



where <5^ ata is again the experimental uncertainty and each [S] is given by Eq. ([8]). The weight 
vector w that minimizes the above \ 2 1S given by w m ; n = — B _1 a, where 

^■=E^ + A ^ (12) 
_ ^ D[ (X k [S ] - X£^) 

D k = x k [sl] - x k [Si } (M) 

The level of agreement between the data set {x data } and the given set of PDFs is now quantified 
- not only by x 2 {X data } — but also by the length of the weight vector w m i n . If |w m i n | < 1 the 
new data set could be included to the original fit within the confidence criterion determined in the 
analysis. That is, the "penalty term" Ax 2 X^fc w k ™ PIP remains below the acceptable value 
Ax 2 - On the other hand, if |w m ; n | > 1, notable tension between the new and the old data is bound 
to exist. 

Applying this method to the neutrino data we use the nuclear PDFs defined in Eq. ([5]) to 
calculate the cross-sections. Therefore, the penalty term splits in two pieces 

15 37 

AX 2 A AEPS09 W k + A *CTEQ6.6 Yl W *> ^ 

k k=l fc=16 
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where Ax^teq = anc ^ EPS09 Ax^p S09 = 50. The key results are given in Table [J for each 
data set separately. The Xw=o ^ s the- value calculated with the central sets, whereas Xw - n * s 
the corresponding value at the minimum of Eq. (jlip . The penalty columns indicate the growths 
induced in EPS09 and CTEQ6.6 and the results are given with and without the normalization 
procedure of Eq. (jlj) . 

The left block of Table U corresponds to the full analysis with all EPS09 and CTEQ6.6 error sets. 
As expected, the normalization improves the x 2 values and diminishes the induced penalties which 
clearly stay within the allowed range. That is, the normalized neutrino data could be included 
in these global fits without an obvious disagreement with the other data. However, without the 
normalization the NuTeV data induces a penalty in EPS09 which starts to get close to the upper 
limit Ax|pg g = 50. Indeed, had we taken the free proton PDFs as fixed (w^ = for k = 16 . . . 37) 
as in [a], the EPS09 penalty would have been much larger. This is demonstrated in the right-hand 
part of Table[IJ Whereas the CHORUS and CDHSW stay well inside the permitted region, the the 
NuTeV data now cause excess penalty in EPS09. That is, there would be a possible contradiction. 

In conclusion, we have demonstrated that disposing the overall normalization by dividing the 
data by the integrated cross-section in each neutrino energy bin separately, all large-Q 2 neutrino 
data show practically identical nuclear effects, consistent with the present nuclear PDFs. Our 
numerical consistency test based on the Hessian method of propagating uncertainties confirms 
that these data could be included in a global fit without causing disagreement with the other data. 

In contrast, without the normalization procedure the nuclear effects preferred by different data 
sets become much more scattered. Especially the NuTeV data seems to display tension with the 
other data. Such is not completely unexpected as in [8| sizable differences in the normalization of 
the NuTeV data among different neutrino energy bins were found. This likely explains the findings 
of f| where, however, all neutrino data were rejected as incompatible. The analysis reported here 
suggests that such a strong conclusion is not justified and we propose a method to deal with the 
apparent tension in different data sets so that the neutrino data can safely be used in global fits. 
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